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Abstract. Modeling the atmospheres of SNe la requires the solution of 
the NLTE radiative transfer equation. We discuss the formulation of the 
radiative transfer equation in the co-moving frame. For characteristic ve- 
locities larger than ~ 2000 km s , the effects of advection on the synthetic 
spectra are non-negligible, and hence should be included in model calcu- 
lations. We show that the time-independent or quasi-static approximation 
is adequate for SNe la near maximum light, as well as for most other as- 
trophysical problems; e.g., hot stars, novae, and other types of supernovae. 
We examine the use of the Sobolev approximation in modeling moving at- 
mospheres and find that the number of overlapping lines in the co-moving 
frame make the approximation suspect in models that predict both lines 
and continua. We briefly discuss the form of the Rosseland mean opacity 
in the co-moving frame, and present a formula that is easy to implement 
in radiation hydrodynamics calculations. 
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There are many astrophysical systems that require a solution to the radia- 
tive transfer equation in moving media; e.g., Wolf-Rayet and other hot stars 
with stellar winds, novae, supernovae, the material surrounding quasars, 
and even the early phases of the universe when the material is still op- 
tically thick (Mihalas 1980). Because there is such a large simplification 
in the radiation-matter interaction terms, it is both customary and expe- 
dient to solve the radiation transfer equation in the co-moving frame. In 
a series of papers Mihalas and co-workers (Mihalas, Kunasz, & Hummer 
1975, 1976a,b,c, Mihalas & Kunasz 1978) examined methods for solving 
the co-moving frame line-transfer problem, where the Doppler effect domi- 
nates because the characteristic width over which the line profile varies is 
small. This effectively increases the importance of the Doppler effect over 
the other 0(v/c) (advection and aberration) effects by the ratio c/v t hcmn 
where ftherm is the thermal velocity corresponding to the intrinsic Doppler 
line- width. Recently, there have been increasingly sophisticated attempts to 
model the atmospheres of hot stars (e.g., Werner 1987), novae (Hauschildt 
et al. 1994, 1995), and supernovae (Branch et al. 1991, Eastman & Pinto 
1993, Hoflich 1995, Nugent et al. 1995, Baron et al. 1995), including NLTE 
for lines and continua, and the effects of line blanketing. Here we system- 
atically discuss the important effects that must be included when solving 
the radiative transfer equation in the co-moving frame, elucidate the range 
of applicability of Eulerian approximations such as the Sobolev approxima- 
tion, and present a co-moving formulation of the Rosseland mean opacity. 
Some of the results discussed here are also discussed in Baron, Hauschildt, 
h Mezzacappa (1995). 



2. Radiative Transfer Equation 



The co-moving frame radiative transfer equation for spherically symmetric 
flows can be written as (cf. Mihalas & Mihalas 1984): 
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We set c = 1; /3 is the velocity; and 7 = (1 — /3 2 ) -1 / 2 is the usual Lorentz 
factor. We emphasize that, in Eq. 1, the physical (dependent) variables are 
all evaluated in the co-moving Lagrangian frame. However, the choice of 
independent variables is free, and the coordinate r in Eq. 1 is an Eulerian 
variable (for a discussion of this point, cf. Mezzacappa & Matzner 1989). 
This is the most convenient choice for solving the transfer equation, where 
one usually specifies the grid by fixing the optical depth for some reference 
frequency. However, this grid differs from the fully Lagrangian grid typically 
used in radiation hydrodynamics. In the latter case, r = r(m). 

In order to illuminate the physics, and without loss of generality, we 
expand Eq. 1 in powers of [3 and keep terms only to O(0). While this is not 
necessary (Mihalas, Kunasz, & Hummer 1976b, Mihalas 1980, Hauschildt 
1992a), it is adequate for most astrophysical flows. To 0(/3), the radiation 
transport equation becomes: 



In writing Eq. 2, we have retained the first term, which accounts for the 
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explicit time dependence of the radiation field in the co-moving frame. We 
have also retained the acceleration term, Both terms are of 0{(3) when 
compared to other terms in the equation, such as the (3/r terms, and hence, 
are of 0{(3 2 ) on a fluid flow timescale and can be dropped (Castor 1972, 
Buchler 1979, Mihalas 1980, Mihalas & Mihalas 1984). Upon doing so, one 
derives the time-independent (or quasi-static) transfer equation in the co- 
moving frame: 



(M + /3)^ + (l-M 2 ) 
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To further simplify the equation and to help elucidate the fundamental 
physics, let us restrict ourselves to consideration of homologous flows: j3 oc r. 
In this case Eq. 3 becomes: 
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In order to identify the physical significance of the terms, it is useful to 
compare this equation to its static counterpart: 
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Comparing Eqs. 4 and 5, the physical meaning of the terms is appar- 
ent: /3^r is the advection term, P/r-Jjj^ represents the Doppler shift, and 
—3(3/rI u describes the effect of aberration. 

It is clear that all three terms are of 0(/3) and must be retained to have 
a consistent treatment in the co-moving frame. The 0(/3) transport equa- 
tion is more difficult to solve because the characteristics, which are simply 
parallel lines of constant impact parameter when the advection term, /3j$jr, 
is neglected, become curved lines; therefore, the reflection symmetry is lost 
(Mihalas, Kunasz, & Hummer 1976b, Mihalas 1980). This is because the 
material is moving, "sweeping up" radiation, causing the characteristics to 
be curved. Therefore, one can no longer use reflection symmetry to inte- 
grate the solution only along outgoing rays. One must integrate along both 
incoming and outgoing rays (Mihalas, Kunasz, & Hummer 1976b, Mihalas 
1980, Hauschildt 1992a). Mihalas, Kunasz, & Hummer (1976b) examined 
the magnitude of the advection and aberration terms and estimated that 
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they are of order 5/3 and that the advection term is more important than 
the aberration term. 

Let us now examine the moments of Eq. 4. The zeroth moment is: 

dJ v 1 d(r 2 H u ) 3 dJ u 

+ a H (3^ — — )=r] u -XuJu, (6) 

or r A or r omv 

and the first moment is: 
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The Eddington moments are given by: 
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where fi = cos 9. When advection is neglected, the moment equations be- 
come: 

ld{r 2 H u ) 3 8J 
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i.e., the gradient of the energy density is absent from the zeroth moment 
equation, and the divergence of the flux is no longer included in the first 
moment equation. Integrating the moment equations (Eqs. 6 and 7) over 
frequency, and assuming that radiative equilibrium holds, i.e., that energy is 
conserved or total emission equals total absorption [J^°(rj v ~XvJv)du = 0], 
one obtains for the zeroth moment: 

Dt r 2 or r 



where we have restored the time derivative: 
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It has been suggested (Eastman & Pinto 1993) that one can correct 
for neglecting the advection term and include the effects of the radiation 
field time dependence on a radiation flow timescale by using Eq. 11 and 
by arbitrarily setting the co-moving luminosity (= r 2 H) to be constant. In 
this case, Eq. 11 becomes: 

Bl = -Mj. (13) 
Dt r y ' 

This is interpreted as an operator equality: 

Dt r v ' 

When Eq. 14 is substituted into the radiation transport equation, the trans- 
port equation becomes: 

ai v , (i-» 2 ) di v (3 r ai v \ 

Comparing Eq. 15 to Eq. 4, we see that this scheme is equivalent to making 
the quasi-static approximation, neglecting advection, and changing the sign 
and coefficient of the aberration term, which is unphysical. 

We have compared simulations with and without the advection term 
for two models that provide reasonable fits to SN 1987A at 13 and 31 
days after explosion. These calculations were performed using version 5.5.9 
of the general radiative transfer code, PHOENIX, developed by Hauschildt 
(Hauschildt 1992a,b, 1993). This code accurately solves the fully relativistic 
transfer equation, Eq. 1, in the quasi-static approximation, = (Note 
contrary to what is incorrectly stated by Hoflich et al., in this volume, we 
do not solve a simple non-relativistic transport equation with the Rybicki 
method, but we indeed solve the full special relativistic, spherically sym- 
metric radiative transfer equation for lines and continua with an operator 
splitting scheme based on a short characteristic method with non-local, 
adjustable approximate A-operator). The model parameters are given in 
Table 1 (for a discussion of the model parameters, cf. Baron et al. 1995). 
For the day-13 spectrum, Figure 1 compares the spectra of a calculation 
that includes advection with one that does not, while Figure 2 displays the 
same for day-31. The differences are about the size predicted by Mihalas, 
Kunasz, & Hummer (1976b), with the effects being more apparent in the 
faster day-13 model than in the much slower day-31 model. Figures 3 and 4 
show comparisons of the temperature profiles for both models. Neglecting 
advection alters the temperature structure, which can be interpreted as 
resulting from the change in the relations between the moments (compare 
Eqs. 6 and 7 with Eqs. 9 and 10). Figure 5 illustrates that this is the most 
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Day 13 Day 31 



T eS (K) 5400 4400 

w (km s- 1 ) 5500 1700 

Ro (cm) 6.1 x 10 14 4.5 x 10 

N 6 6 



TABLE 1. T off , Ro , v are the "effec- 
tive" temperature, radius, and velocity 
at the reference point r st d = 1; N is the 
power-law density index. 



important effect of neglecting advection. The temperature structure and 
departure coefficients are kept fixed in order to only alter the transport 
equation. In this case, the emergent spectra are much more similar than 
those in Fig. 1. From these results, it is clear that the effects of advection 
should not be neglected in models of supernovae where the characteristic 
velocities are larger than ~ 5000 km s _1 . For systems where velocities are 
lower than ~ 2000 km s _1 , advection may be neglected with reasonable 
accuracy. 



3. Quasi-Static Approximation 

We can estimate the effects of making the quasi-static approximation by 
examining the radiation transfer equation. For clarity, we again restrict 
ourselves to O(P) and homologous flows. Restoring the time derivative in 
Eq. 4, we have: 

d 4 + o. + «^ + "-"'>"- 



dt dr r dfi 



r 1 d In v 

= Vv- Xuh- (16) 



In order to solve this equation numerically, we would replace the time 
derivative with the difference: 

f)T T — T n 

dt ~ st ' { ' 

where I v n is the intensity evaluated at the previous time t n , I u is the inten- 
sity at the current time t n+1 , and 5t = t n+l — t n . Inserting this expression 
into Eq. 16, and moving the time derivative to the right-hand side of the 
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Figure 1. The spectrum produced by a full calculation, which fit SN 1987A at 13 
days past explosion (solid line), is compared to one with the same parameters in which 
advection is neglected (dot-dashed line). Both calculations are in radiative equilibrium. 



transfer equation, we obtain: 

or r Ofi r { o in v 

j n i 

= ^ + + (18) 

which shows that the time derivative term can be viewed as an additional 
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Figure 2. The spectrum produced by a full calculation, which fit SN 1987A at 31 
days past explosion (solid line), is compared to one with the same parameters in which 
advection is neglected (dot-dashed line). Both calculations are in radiative equilibrium. 

source and sink of radiation. We can estimate the size of the error made 
in the quasi-static approximation by examining the ratio x _1 /(c<5i), where 
we have restored the explicit c. In supernovae, the natural timescale is the 
age of the object, t = R/v. We may estimate that x ~ T /R, where r is an 
appropriate optical depth. Then, the ratio becomes: 



(19) 
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Figure 3. The temperature profiles for the models plotted in Fig. 1 are compared. 



For Type la supernovae at maximum light, the continuum extinction opti- 
cal depth is about 10, and (3 ~ 1/30. So the error is at most 0.3%; small 
compared to errors in the atomic physics. This error will be considerably 
smaller for other types of supernovae, which are optically thick for longer 
times. In fact, Eq. 19 is an overestimate of the error because we have ne- 
glected the source term I v n /(c5t), which counteracts the extra sink term. 

Claims (Eastman, this volume) that in order to account for "old pho- 
tons" , time dependence must be included in the transfer equation are not 
correct. The effects of both dynamic and static diffusion are included in ra- 
diation hydrodynamics without including the time-dependent term in the 
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Figure 4- The temperature profiles for the models plotted in Fig. 2 are compared. 

transfer equation (Mihalas & Mihalas 1984). The effects of departures from 
radiative equilibrium can be included in our modeling. 

4. NLTE Effects 

Figure 6 displays the model atoms for Li I, Ca II, Ti I, Ti II, Fe II, and 
Co II used in PHOENIX. In addition H I, He I, He II, Mg II, Ne I, and O I are 
also treated in NLTE. In the very near future we will add C I-IV, N I- VI, 
O II-VI, Si II— III, and S II— III. In particular, our large Fe II model atom 
(617 levels, 13675 primary transitions; see Hauschildt & Baron 1995, for 
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Figure 5. The spectrum produced by a full calculation, which fit SN 1987A at 13 days 
past explosion (solid line) , is compared to one with the same parameters in which advec- 
tion is neglected (dotted line). The structure is fixed to be that of the full calculation. 

more details) allows us to determine the importance of the size of model 
atoms and of the validity of common assumptions that are often made in 
handling the millions of secondary transitions that must be included in 
order to correctly reproduce the UV line blanketing. We find that large 
changes in luminosity (see Pinto, this volume) can be avoided by treating 
the secondary lines with a small but non-zero thermalization parameter 
(Baron et al. 1996), as required by physical considerations. 




NLTE Modeling of SNe la Near Maximum Light 13 




Figure 6. Grotrian diagrams for Li I, Ca II, Ti I, Ti II, Fe II, and Co II model atoms used 
in the calculations. Only the "primary" radiative transitions are displayed; see Hauschildt 
& Baron (1995) and Hauschildt et al. (1996) for more details. 
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5. Sobolev Approximation 

The Sobolev approximation developed by Sobolev (1960) and Castor (1970) 
and extended by Hummer & Rybicki (1985, 1992) and Jeffery (1989, 1990, 
1995, 1996) has proven extremely valuable in providing line identifications 
and minimum and maximum velocities in supernovae (Branch et al. 1983, 
1985, Jeffery & Branch 1990, Jeffery et al. 1991, 1992, 1994, Kirshner et al. 
1993, Filippenko et al. 1992), because it allows one to calculate line profiles 
without solving the transfer equation, it is nearly analytic, and quite conve- 
nient. The above analyses were concerned with identifying strong lines, and 
continuum effects were neglected. More recently the Sobolev approximation 
has been used to solve the rate equations for detailed model atoms including 
continua (Eastman & Pinto 1993, Hoflich 1995, Hoflich et al., this volume). 
However, because the escape probability is derived by neglecting the effects 
of neighboring lines, it is only valid for isolated lines, and is invalid when 
there are many weak overlapping lines (Castor 1970, Rybicki 1984, Avrett 
& Loeser 1987). This is likely to be the case in the UV, where line blan- 
keting is severe. Rybicki (1984) also has discussed that escape probability 
methods such as the Sobolev approximation are inaccurate at small line 
optical depths, particularly when there are many overlapping lines. Since 
the source function predicted by the Sobolev approximation at the surface 
is incorrect by a factor of y/e, where e is the line thermalization parameter, 
the value of J found by the formal solution will also be in error, which in 
turn will lead to errors in the rate equations. In Figure 7 we display the 
number of overlapping lines in a range of 6 intrinsic Doppler widths around 
any given wavelength, as a function of wavelength, for the day-13 SN 1987A 
model, which has a statistical or micro-turbulent velocity of £ = 50 km s -1 . 
The lines are said to overlap if, for any particular line, another line has its 
line center ±6 intrinsic line-widths from the reference line. The Doppler 
widths are calculated at deepest depth point in the model. As expected, in 
the UV the mean number of overlapping lines is typically around 100, and 
can be as large as 500. This implies that the radiative transfer in SN (and 
nova) atmospheres must explicitly include the effects of overlapping lines 
and continua. Otherwise, the radiative rates for these transitions would be 
incorrect, particularly in the outer parts of the atmosphere where ionization 
corrections are most important. Although this requires a very fine wave- 
length grid for the model calculations, detailed models can be computed 
using modern numerical techniques on even moderately sized workstations. 
Thus, the Sobolev approximation cannot be used in detailed NLTE calcu- 
lations for SNe, because the radiative rates calculated in this approxima- 
tion are inaccurate. Similar results have been obtained by Hauschildt et al. 
(1995) in nova model atmosphere calculations. 
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Figure 7. The number of overlapping lines in a range of 6 Doppler widths around any 
given wavelength, as a function of wavelength, for the day- 13 SN 1987A model. Only 
lines stronger than 10~ 3 times the local b-f continuum are included. Clearly, the Sobolev 
assumption that individual lines do not overlap is not fulfilled. 

For pure hydrogen atmospheres, Duschinger et al. (1995) found good 
agreement between the non-relativistic Sobolev approximation and non- 
relativistic co- moving frame full transport calculations, which shows that 
the Sobolev approximation is accurate for well separated lines such as the 
Balmer lines. However, this situation is not reproduced in most spectral 
regions, and therefore, the Sobolev approximation is of limited use for de- 
tailed modeling of SN or nova envelopes. 

6. Expansion Opacities 

Radiation hydrodynamic calculations of supernova light curves require ac- 
curate fluxes, and it has long been realized that the static Rosseland mean 
opacity does not produce an accurate flux in moving atmospheres (Karp 
et al. 1977). The work of Karp et al. (1977) provided an approximate for- 
mula for the Rosseland mean opacity in the observer's frame. However, 
nearly all radiation hydrodynamics calculations are performed in the co- 
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moving frame; hence, a co-moving formulation is required. 

We have derived the Rosseland mean opacity in the co-moving frame to 
0(/3) (Baron, Hauschildt, & Mezzacappa 1995). Let us first recall that the 
static Rosseland mean, x% is given by: 

Xr k Jo dT 

To derive the non-static Rosseland mean, we will assume homologous flows. 
In addition, we make the Eddington approximation, K u = 1/3 J^, implying 
that the co-moving radiation field is close to isotropic, which is an excellent 
approximation in the diffusive regime (large optical depth) because the 
radiation is collision dominated (Pomraning 1982). 

We find that the co- moving Rosseland mean opacity to 0([3) is given 

by: 

^-(x+^r 1 (2i) 



in the gray case, and: 



Xr 



Xr 



XrT r 7r jo ol a In v 



(23) 



in the non-gray case. In deriving Eq. 23, we have used ^ = 4(a / tt)T 3 ^ . 

It follows that the co-moving multi-group flux to be used in radiation 
hydrodynamics is given by a Fick's law diffusion equation: 

We emphasize that \v m Eq. 20 and 23 contains contributions from con- 
tinua, lines, and scattering opacities, and nowhere have we had to treat 
lines differently from continua. 

In the case that the opacity may be approximated by a power-law, 
Xu oc v~ n , the last integral in Eq. 23 may be evaluated by an integration 
by parts, yielding: 

r it Jo oToniv vxr 
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and therefore: 

^^-(« + 2)V (26) 

We have calculated the correction factor for atmospheres appropriate to 
Type II supernovae. As an illustrative case, Figure 8 displays the density, 
temperature, x% Xr-> an d the effective value of n (that is, the value of n 
one obtains from Eq. 26 using the exact values of x°r an d Xr) as functions 
of r for the day-13 model. As expected from Eq. 21, the largest correction 
occurs at low optical depth (the formula breaks down at very small optical 
depths since the assumptions used to derive it, i.e., LTE and the Eddington 
approximation, are not fulfilled), and the correction is essentially irrelevant 
at high optical depths, where 1/x r << ft/r. 

7. Conclusions 

We have shown that advection cannot be neglected in the co-moving so- 
lution of the radiation transport equation. Its main influence is on the 
temperature structure, through the term it adds to the equation of radia- 
tive equilibrium. The errors made in neglecting advection scale with the 
velocity; while it may be acceptable to neglect advection in systems where 
the velocities are < 2000 km s _1 , such as in hot stars, novae with low wind 
velocities (e.g, Nova Cas 1993), and Type II supernovae at late times, it 
cannot be neglected for supernovae at early times and novae with high wind 
velocities (e.g., Nova Cygni 1992). 

We have also shown that the Sobolev approximation is likely to be 
invalid for weak lines in the co-moving frame, since many of these lines 
overlap. 

We have derived an approximate expression [good to 0(/3)] for the 
Rosseland mean opacity that can be used in radiation hydrodynamics cal- 
culations. The Doppler shift is fully accounted for in this approximation. 
Our formula shows that, at large optical depths, the static Rosseland mean 
is accurate, and hence, for all radiation hydrodynamics calculations that 
use flux-limited diffusion, the static approximation is excellent. 
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Figure 8. The effective value of n, the density, the temperature, Xr (solid line) , and \ R 
(dashed line), as functions of r st d, for the day-13 model. At very low optical depths, the 
formula breaks down. 
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